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^ ■ Abstract 



We propose a novel and simple method to compute the partition function of statistical me- 
chanics of local and semi-local BPS vortices in the Abelian-Higgs model and its non-Abelian 
extension on a torus. We use a D-brane realization of the vortices and T-duality relation to 
domain walls. We there use a special limit where domain walls reduce to gas of hard (soft) 
one-dimensional rods for Abelian (non-Abelian) cases. In the simpler cases of the Abelian-Higgs 
model on a torus, our results agree with exact results which are geometrically derived by an 
explicit integration over the moduli space of vortices. The equation of state for U(N) gauge the- 
ory deviates from van der Waals one, and the second virial coefficient is proportional to 1/y/N, 
implying that non-Abelian vortices are "softer" than Abelian vortices. Vortices on a sphere are 
also briefly discussed. 
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1 Introduction 

Knowledge of moduli space structure of solitons is important in order to understand not only their 
own dynamics but also non-perturbative effects in field theory. The moduli space structure of the 
solitons signifies its topology, metrics and singularities. The volume of the moduli space also pos- 
sesses an important meaning and plays an essential role in understanding the non-perturbative 
dynamics. Recently Nekrasov has shown that the prepotential of M=2 supersymmetric gauge 
theory in d = 3 + 1 can be obtained from some statistical partition function, whose free energy 
measures the volume of the moduli space of Yang- Mills instantons pQ. The instanton moduli 
space is non-compact and of course the volume diverges, but a coefficient of the divergent part 
gives perturbative and non-perturbative information on supersymmetric gauge theory. It is also 
interesting that the partition function or the volume of the instantons is closely related to topo- 
logical string amplitudes on suitable Calabi-Yau manifolds. Precisely speaking, the instantons, 
which are used for the calculation of the volume, are not solutions of the equations of motion 
in supersymmetric gauge theory. However, in the calculation of the partition function or the 
volume, the so-called localization theorem says that the result does not depend on details of 
the moduli space structure. This suggests that we do not need to know the exact solutions and 
metrics in order to calculate the volume of the moduli space of some class of BPS solitons. How- 
ever, the issue has not been settled except for Yang-Mills instantons, since other examples and 
applications have not been investigated yet. For instance, similar method should be applicable 
to calculate the effective (super)potential of d — 1 + 1 supersymmetric gauge theory from the 
statistical partition function associated with the volume of the vortex moduli space. 

Another direct application of the volume of the moduli space is a classical statistical mechanics 
of the solitons. The Gibbs partition function of solitons at finite temperature T is given by an 
integration over a phase space, which is a cotangent bundle T*Ai of the moduli space A4 of 
solitons: 



Here D represents the dimensions of the moduli space, and H(x,p) is a Hamiltonian of the 
soliton system in terms of the moduli parameters. If we assume that the solitons are sufficiently 
diluted and interactions can be ignored, we can regard the Hamiltonian as quadratic in momenta 



the integration over the momenta (cotangent direction). The partition function becomes 



which is proportional to the volume of the moduli space [2]. Many applications and calculations 
from this point of view can be found [2] [6] in the case of Abrikosov-Nielsen-Olesen (ANO) 
vortices [7j, namely vortices in the Abelian-Higgs model at the critical coupling (the BPS limit). 






metric on A4. Then we can explicitly perform 
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On the other hand, non-Abelian BPS vortices have been recently found [El [9] and extensively 
discussed in the (supersymmetric) non-Abelian Higgs model, which is the U(N C ) gauge theory 
coupled with N F (> N c ) fundamental Higgs fields. Like (non-commutative) instantons, single 
vortex in the Nq = N-p case is found to carry internal moduli CP^ 0-1 in addition to position 
moduli. They can confine monopoles in the Higgs phase pm E] an d this fact is applied [T2l [T3] 
to show coincidence of the BPS spectra of <i = 3 + 1, M = 2 supersymmetric gauge theory and 
d = l + l,A/" = 2 supersymmetric CP" model. The moduli space of multiple (k-) vortices has 
been conjectured in string theory [8j and then has been completely determined in field theory 
[H I fl5 l [T6]. which has turned out to be some resolution of fc-symmetric product of C x CP^ " 1 
|14j . Although its explicit metric is still unknown except for a single (k = 1) vortex, an integration 
formula of the Kahler potential has been found [17]. The structure of the moduli space of two 
{k = 2) vortices has been worked out [TH [16] and is applied to non-Abelian duality [19], and 
classical dynamics of non-Abelian vortex-strings has been clarified [20]. The extension to a 
superconformal field theory [21 J as well as the Chern-Simons-Higgs theory [22] has been discussed. 

In this article, we propose a novel and simple derivation of the volume of the moduli space of 
Abelian as well as non-Abelian (semi-)local vortices in order to apply to the statistical mechanics 
of vortex gas. We utilize a geometric and stringy (D-brane) interpretation of the Abelian and 
non-Abelian BPS vortices in supersymmetric gauge theory [8] and use a T-dual relation between 
the vortices and domain walls, which was proposed in [231 El] • This interpretation shows us a 
schematic structure of the vortex moduli space without any detailed information on the exact so- 
lutions and metrics. (See Fig. 17 of Ref. [21] for the moduli space of two non-Abelian local vortices 
in U{2) gauge theory.) It is interesting to observe that the T-duality reduces the calculation of 
the volume into a simple problem of a gas of hard rods between a one-dimensional interval in the 
Abelian case. The volume of the moduli space of the ANO vortices in the Abelian-Higgs model 
(with iV c = N F = 1) has been already calculated by [2] [1] (see also §7 in [5]) and an analogy 
with the hard rod problem was pointed out there. However this relation has been mysterious for 
a long time. We can explain from a point of view of the T-duality in superstring theory why the 
one-dimensional hard rod problem appears. The advantage of our method is that it can be easily 
extended to the general cases, namely local and semi-local, Abelian and non-Abelian vortices. 
We find that the T-duality enables us to reduce the problem of a gas of non-Abelian local vortices 
to a gas of soft rods, in contrast to the hard rods in the case of Abelian local vortices. We obtain 
that the second virial coefficient is proportional to 1/yN for non-Abelian local vortices in U(N) 
gauge theories at large N. This shows that the exclusion volume of a vortex behaves as l/\/N 
as we increases N, getting closer to an ideal gas. Moreover, equation of state for local vortices in 
U(N) gauge theory deviates from van der Waals one, contrary to the case of the Abelian-Higgs 
model. 

Comparing with the calculation in [2]-[5], our result in the Abelian case gives the precise 
answer even though we take a special limit to the configuration and do not use the exact metrics. 
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So we expect that this is an another example where the localization theorem effectively works. 
In fact calculation of the volume of vortex moduli space using the localization theorem can be 
found in the case of the Abelian-Higgs model [6j. 

In Sec. 2 we give a brief review of brane configurations of vortices and domain walls. In sec. 3, 
we calculate the partition function of vortices in the Abelian-Higgs model, which agrees with the 
previous result. In sec. 4, we give general formula of the partition function of local/semi-local 
non-Abelian vortices which is new and shows the power of our method. We perform explicit 
calculation in several cases. Sec. 5 is devoted to a discussion. An interesting duality is observed 
and our method is extended to the case of a base manifold S 2 . Appendix gives some details of 
derivation of the virial coefficient. 



2 Non-Abelian Vortices and D-brane Interpretation 



We first start with (2+l)-dimensional U(Nq) gauge theory with Np(> Nq) massless Higgs fields 
in the fundamental representation. The Lagrangian is given by 



C = Tr 



-^- 2 f mn f mn + v M H {v M H) j - L ( M t _ clNc y 



(2-1) 



where M,N are the indices of space-time (M,N = 0,1,2) and the space-time metric is cho- 
sen as T]mn = diag(+l, — 1, — 1). The field strength is defined by Fmn = — i\PM, T^n\ = 
d M W N - d N W M + i[W M , W N ), where W M is the U(N C ) gauge field. The scaler fields H rA (r = 
1, • ■ • , Nc), (A = 1, • • • , iVp) are expressed by elements of an Nq x N-p matrix and the covariant 
derivative is defined by T> M H = OmH + iWmH. The constants g and c (> 0) are the gauge 
coupling constant and the Fayet-Iliopoulos (FI) parameter respectively. 

For static configurations, the energy is bounded from below as follows: 



E 



d xTr 



j 2 (f 12 + ^(c1 Nc 



HH*)\ +V S H(V S H? -cF 12 



> -c / d 2 xF, 



12 



luck, 



(2.2) 



where the integer k = / d x F 12 represents the number of vortices (vorticity). Since the 

2vr J 

configurations of the BPS vortices saturate the inequality Eq. (12.21) . they satisfy the following 
BPS equations (the vortex equations): 



V S H = 0, 



F 12 + ^-(cl Nc 



HW 



0. 



(2.3) 



If the fluctuation energy above the energy of the multi- vortex static configuration is less than that 
of the mass gap between massless moduli fluctuations and massive fluctuations, the dynamics of 
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the system is well described by the geodesic motion on the moduli space Aik of the k vortices 
(Manton's moduli/geodesic approximation [25]), so that the Lagrangian for multi- vortex system 
is given by 

L = g i3 (<l>,4>)<j> i fi, (2.4) 

where (f) 1 , <^*, ( % = 1, ■ ■ ■ ,D = dime A4k) are complex coordinates of the moduli space and 
(?jj(0,0) is a Kahler metric of M.}.. 

In the cases of compact base manifolds, there is the maximum number for vortices allowed 
to exist [26]. For general Nq with a torus T 2 , this maximum number can be obtained by taking 
the trace of the second equation of the BPS equation ( 12. 3ft and by integrating over T 2 as 

2 2 r 

-2nk + N c ^A = 9 — I d 2 xTr (M f ) > A > jfeBjvb, (2.5) 

where A is the area of the torus and we have defined the Bradlow area 

1 4vr . , 

B "c = TT — • 2.6 

JMc g c 

The limit to saturate this inequality is called the Bradlow limit. In particular, we can regard 
B\ = ^p- c as an effective area of the Abelian (ANO) vortex. Roughly speaking, the inequality 
( 12.51) implies that we can squeeze as many vortices as Nc times the maximal number of the 
Abelian (Nc = 1) vortices for a given area A of the torus. This behavior can be understood 
intuitively as due to the freedom for vortices to avoid occupying the same points in their additional 
internal moduli space CP^ 0-1 when they are squeezed too much in the actual configuration 
space. This freedom allows the non- Abelian vortices to overlap in the configuration space more 
easily compared to Abelian vortices. Therefore the effective area of the non-Abelian vortex B^ c 
becomes 1/iVc < 1 times the area of the Abelian vortices B\. The Bradlow area Bn c indicates 
the smallness of the effective area of vortices at the limiting high density of vortices. Note that 
there exists another inequality in the case of non-Abelian gauge theories with iVc > k. Since 
all vortex solutions in this case can be embedded to a theory with a smaller number of colors 
Nq(= k') = k, the vortex configuration for N c > k should satisfy the inequality for Nq = k, that 
is, 

A>^ (=3,) (2.7) 
9 c 

stating that the area of the torus must be larger than that of one ANO vortex. 

As we mentioned in Eq. (11.21) . the classical partition function of k- vortex system is propor- 
tional to the volume of the moduli space Aik- To calculate this partition function (volume), 
we utilize the duality relation between vortices and domain walls [23l [24] . If we compactify re- 
direction on S 1 as x 2 ~ x 2 + 2imR, a vortex configuration can be viewed as domain walls. The 
profile of the kink solution of the domain walls is described by the eigenvalues of S defined by 



^ 1) = -2^ l0S 



Pexp^i^ dx 2 W 2 {x 1 ,x 2 )\ 



"2irR 

{21 
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(a) (b) (c) 

Fig. 1: Duality between vortices and domain walls: Nq = Np = 1, k = 2. (a) Profile of energy density 
of 2-vortex configuration on IxS 1 with period 2ttR, (b) Profile of ^{x 1 ) = f dx 2 W2, (c) Profile 

of ^(x 1 ) written on R x S 1 with period 1/R. 



where P stands for the path-ordered product. Note that under the gauge transformation U = 
e - mx2 / R 1 Nc (ji g g), E^ 1 ) transforms as E — > E + n/Rl^ c and there is an identification 
E ~ E + n/RX]sr c . Thus the eigenvalues of ^(x 1 ) are interpreted as a function which takes value 
in S 1 with radius 1/R. Fig. [1] depicts an example of the Nq = Np = 1 case. The k- vortex 
configuration (Fig. [H-(a)) corresponds to the /c-wall configuration (Fig. [H-(b)) in the fundamental 
region, and it also can be interpreted as the configuration of E^ 1 ) with k windings (Fig. [U-(c)) 
on the cylinder. 

This duality can be regarded as a T-duality between the brane configurations of vortices 
[8] and domain walls [28]. Our model can be embedded into the supersymmetric system with 
eight supercharges and the associated model is realized by a combination of various kinds of 
branes in Type IIB superstring theory. The brane configurations are expressed in Table 1 and 
drawn in Fig. 2 schematically. The (2 + l)-dimensional U(Nq) gauge theory coupled with N-p 
massless hypermultiplets is realized on iVc D3 branes in the Hanany-Witten setup as in 
Fig. [2- (a), since the x 3 direction of the D3 brane worldvolume is a finite line segment. Dl 
branes correspond to vortices in the D3 brane worldvolume theory, since they are interpreted as 
co dimension- two objects on the D3 branes [8]. Fig. [2] shows the brane configuration T-dualized 
along the x 2 -direction. The theory without vortices is mapped to the world-volume theory on 
the Nq D2 branes, which is the (1+1) dimensional U(Nc) gauge theory with N-p (massive) 
hypermultiplets. The vortices are mapped to the kinky D2-branes representing domain walls 
[271 [28]. The eigenvalues of E^ 1 ) can be interpreted as the position of D2 branes in S 1 . 
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(a) Brane configuration for vortices 
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(b) Brane configuration for domain walls 



Table 1: Brane configurations for vortices and domain walls: Branes are extended along directions 
denoted by •, and are not extended along directions denoted by — . The symbol x denotes the codi- 
mensions of the k Dl-branes (D2'-brane) on the worldvolume of the D3-branes (D2-branes) excluding 
the x 3 -direction which is a finite line segment. 

3 A Limit of the Profile and the Moduli Integration: the 
Abeain-Higgs Model 

If we use the T-dual relation via the brane configuration, the evaluation of the volume of the 
vortex moduli space reduces to a calculation of the volume of the domain wall configurations 
(kink profiles). However, it is difficult to solve the domain wall equations and to integrate over 
their configuration space in general. 

In order to proceed the evaluation, we demand an approximation which simplifies the profile 
of the domain wall solution. Now, let us consider a special limit 

2 

g 2 c, 1/R— > oo, with d = — — — : fixed. (3.1) 

g 2 cR 

Of course, this is a rough approximation with respect to the kink profile, we will amazingly find 
that we can obtain the exact results even in this limit. Note that using d defined above the 
inequalities (12.51) and (12. 7p are translated in the T-dualized picture to 

2ttR' > —kd , 2txR' > d, (3.2) 



6 



X 



1,2 



NS5 NfD5 ^ 



NS5 



4,5, 6 




NcD3 





















£ 


Nc 


D3 



Dl 



(a) Brane configuration for vortices 




(b) Brane configuration for domain wall 



Fig. 2: Brane configurations for vortices and domain wall 



respectively. These inequalities can be easily shown even in this T-dualized picture as dicussed 
below. This fact is the first evidence for the usefulness of our T-dualized picture in this paper. 

Here we describe a limit shape of the ANO vortices in the Abelian-Higgs model (Nq = Np = 1) 
for simplicity. In this case, the moduli space of the multi-vortex system consists of the position of 
each vortex. The positions of vortices in the x 2 -direction are translated into the phase degrees of 
freedom of domain walls after the T-duality along x 2 direction. If one vortex is located at z — 0, 
the corresponding function ^(x 1 ) is given by 



tlx 1 ) = I 4- 



i 

Rd 





, d 

X+ 2 
1 

R 



for 
for 
for 



x 1 < 
d 



d 
2 



— < x 1 < - 
2 2 



(3.3) 



x 1 > 



d 



From this expression, we find that the parameter d represents the effective thickness of the dual 
domain wall. Fig. [3] shows the profiles of ^(x 1 ) before (Fig. [3]-(a)) and after (Fig. (3]-(b)) taking 
the limit. We have previously shown both from field theory [30] and string theory (the s-rule) 
[28] that these domain walls cannot pass through nor be compressed with each other. Therefore 
these objects can be regarded as 1-codimension rigid bodies with length d. Using this fact, we 
can map the multi-vortex system to the gas of hard rods, namely 1-dimensional gas of the rigid 
bodies. 

Using the above limit, where the domain walls are regarded as the hard rods in 1-dimension, 
let us now consider the vortices on a torus T 2 with periods 2ttR and 2irR' and dual domain wall 
configuration space. Fig. HJ-(a) shows the profile of ^(x 1 ) and Fig. 0]-(b) shows the corresponding 
system of hard rods. For the gas of identical hard rods with mass m on S* 1 with period L, the 
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(a) (b) 
Fig. 3: Profiles of ^(x 1 ) before and after taking the limit g 2 c, 1/R — ► oo with d = 2/g 2 cR finite. 
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(b) 1-d gas of hard rods on S 1 



Fig. 4: Profile of ^(x 1 ) for vortices and corresponding 1-d gas of hard rods. 



classical partition function can be easily calculated as 

k 



J rods 



1 

jfe! 



(3.4) 



In this case, the period is L = 2irR! and each rod has mass m = 2nc which corresponds to 
the mass of the vortex. There are additional phase degrees of freedom which correspond to the 



positions of vortices in x 2 -direction. In the limit g 2 c, 1/R — > oo with d 



g 2 cR 



finite, these phase 



degrees of freedom become independent and have the period 2nR. In other words, the T-duality 
maps the Fl-parameter c to 2ttRc. So we should replace c by 2-kRc in the vortex picture. (Note 
that the combination d = is invariant under the T-duality.) Therefore the partition function 
of fc-vortex system on a torus T 2 is given by 



yN C =N F = l 

^k,T 2 



i (cT) k 2ttR' (2nR' - kdf- 1 (2nR) k = 1 (cT) k a(a~^ 
k\ k\ \ g z c 



k-l 



(3.5) 



where A = (2n) 2 RR' is the area of T 2 . This result coincides with the exact partition function 
which has already been known HI IS]- The inequality (12. 5p with iVc = 1 implies that the 
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number of vortices k must be less than A In the Bradlow limit k = ^A, the partition function 
Eq. (13. 5p vanishes. In our interpretation of vortices as 1- dimensional rods, this maximum number 
can be understood from the fact that sum of the length of rods can not exceed the compact period 
of S 1 , namely -f^ = kd < 2ttR' from the first equation of (I3.2p . 

We can derive the van der Waals equation of state for the vortex gas in the thermodynamic 
limit A — > oo, k — > oo with keeping k/A fixed (see eg. §7.15 in [5]), 

Here P is the pressure of the vortex gas and the Boltzman constant is unity (ks = 1) in our 
notation. From this we find the pressure of the vortex gas diverges at the maximal number 
of vortices determined by the area and couplings. Compared with the general van der Waals 
equation of state (P + ak 2 /A 2 ) (A — kb) = kT we have a = and b = Arc/ {g 2 c). The former can 
be understood as the BPS property that there exists no potential energy between BPS vortices. 
The latter implies the size (exclusion area) of the ANO vortex is 4ir/(g 2 c). 



4 Local/Semi- local Non-Abelian Vortices 

So far we have seen that the partition function for multi-vortex system can be calculated by 
using a system of hard rods in the case of the Abelian-Higgs model (iVc = Np = 1). We next 
extend this method to the more general cases of non-Abelian and Abelian gauge theories: local 
non-Abelian vortices (Nq = Np > 1) and semi-local (non-Abelian) vortices (Np > Nq). These 
cases were not known previously, which proves the power of our method. To this end, we treat 
vortices in a model with a twisted boundary condition 

H(x\ x 2 + 2ttR) = H(x\ x 2 )e 2mRM , (4.1) 

where M = diag(«i/i?, • • • , ajv P /i2), a\ < ■ ■ ■ < «at p < 1 + &x- We can reproduce the partition 
function with ordinary boundary condition by taking the limit — > after calculating the 
partition function with the twisted boundary condition. In this case, fc-vortex configuration 
corresponds to kN-p domain walls represented as kinks of the eigenvalues of ^(x 1 ) as shown in 
Fig. [5l By taking the limit (13.11) . we can identify each domain wall with a 1-dimensional rod as 
before. In this limit, all kinks of ^(a; 1 ) which correspond to domain walls have the same slope, 
which is 1/dR = g 2 c/2. In this case there are Np types of rods whose masses are given by 

m A = 2Trc(a A +i - a a), A = 1, ■ ■ ■ , N F - 1, 
m Np = 2nc(l + oc\ - a Np ). (4.2) 

and the period of phases are given by 

I A = 2nR/(a A+1 - a A ), A = 1, • • • , N F - 1, 
l Np = 2nR/(l + ai -a NF ). (4.3) 
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Fig. 5: Eigenvalues of ^(x 1 ) as the dualized configuration of k vortices. Kinks represent kNp domain 
walls. 



Then we can calculate the partition function for the multi-vortex system in terms of the 1- 
dimensional hard rods. Note that we have to take into account the fact that there are k indis- 
tiguishable sets of rods corresponding to k indistinguishable vortices. Therefore the partition 
function should be divided by k after integrating over the configuration space of distinguishable 
rods. Then the partition function for multi-vortex system takes the form of 

, rp n kN F / N ¥ 

Z?°f* = l^—J \T[mAlAj -j- J^dx 1 ---dx kNF 

—J (2^)^(2^)^- J d Vl ... dy kNF . (4.4) 

Here the integration with respect to positions {x n } is taken over the configuration space of the 
rods, which is denoted as X. In the last line, we have redefined the coordinates as y n = x n — oiAd 
if the n-th rod sits between the A-th flavor and the (A + l)-th flavor, and the corresponding 
domain of integration is denoted as X' . The redefinition of the positions corresponds to consider 
a slant torus (Fig. [6]) where domain wall configurations are right-angled and perpendicular to 
the x 1 -direction. Therefore we can associate the configurations of vortices with that of rods with 
length d and particles with zero size pierced by the rods. Recalling the original profile of the kink 
solution, we can see there exist some overlapping rules of the rods; a part of rods which is divide 
by the particles can overlap with each other, but not be allowed for some of these. So we also 
introduce colored rods in order to clarify the exclusion rules of rods. (See Fig. [71) Note that the 
above formula in the last line is completely independent of the parameters {a^}. Therefore, our 
result with the twisted boundary condition (14.11) (non vanishing a a) is applicable to the ordinary 
case with non- twisted boundary condition (at a = 0). 

The partition function Eq. (14.41) can be interpreted as the asymptotic form of partition func- 
tion of vortex gas on the rectangular torus in the limit Eq. (13. ip . Let us assume that the partition 
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Fig. 6: The profile of eigenvalues of ^(x 1 ) and slant configuration for Nq = Np = 2, k = 2. 
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Fig. 7: (a) Configuration of rods ( ■ — ■ and particles (•) pierced by the rods for = Np = 2, k = 2 

corresponding Fig. [H (b) Those for A^c = Np = 3, k = 2. Each rod is divided by particles into the iVc 
parts with different colors. Each part of rods can overlap with parts of the other rods with different 
colors but not with parts with the same color. 

function is independent of details of the torus and depends only on the area of torus A. This is 
the case for the partition function in the model with Nq = Np = 1. Then we can calculate the 
partition function for vortex gas on a surface which is topologically a rectungular torus with area 
A. In addition, we can show that Eq. (14 .4p gives the exact form of partition function as follows: 
The limit Eq. (13.11) can be rewritten as 

From the explicit expression of the metric of moduli space [15] and dimensional analysis, we can 
show that the partition function takes the form of 



fi — » oo, R , c , g : fixed. 



(4.5) 



yN C ,N F 



(cT) kNp A kNp f (g 2 cA) 
(c T) kN -A kN -f(g 2 c A ) 



A = (2tt) 2 R R 



(4.6) 



where f(g 2 cA) is an unknown funtion of dimensionless parameter g 2 cA. From this expression, 
we find that the partition function is independent of \i. Therefore the limit fi —>■ oo gives the 
exact partition function. In the discussion below, we assume that the partition function depends 
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only on the area of torus A. We will check that the partition function Eq. ( 14.4ft indeed gives the 
exact result for one vortex with Nq = Np = N and Np > Nq = 1. 

We now perform the integration in (14.41) explicitly in some cases, all of which are new results. 



1) Two (k = 2) local non-Abelian vortices with iVc = Np = 2 

Let us consider the domain of integration for Nq = Np = 2 and k = 2 for example. The 
domain can be easily seen from Fig. [6] and Fig. [TJ-(a). In this case there are two rods with one 
particle inside each of them. These rods can overlap with each other contrary to the case of 
iVc = Np = 1. However the edges of the rods cannot overlap with the particles inside the other 
rods. Therefore the domain of integration is given by 



X' = < (2/1,2/2,1/3,2/4) 



< 2/1 < 2-kR', 2/1 < 2/2 < 2/1 + d, 2/2 < 2/3, 2/i + d < 2/4, 
2/3 < 2/4 < 2/3 + d, 2/3 + d < 2/2 + 2nR', 2/4 < 2/i + 2irR' , 



(4.7) 



Performing the integral in Eq. ( 14. 4ft over the domain ( 14.7ft . we immediately obtain the partition 
function of the vortices (the volume of the domain wall (rods with particles) configuration space) 



1, rr^A ( '4VT 



2 



(cT)M^f) A[A-~\ for ^<A 

7 n c =2,n f =2 _ I 2 \g 2 cJ V 3(? 2 c/ g 2 c (AR) 

-(cT) 4 I A —J A I — A j for <A< — 

\ 6 V 9 cj \9 c J g 2 c g 2 c 

Note that there exists a lower bound for the area A > An/ {g 2 c) which corresponds to the Bradlow 
limit (12.51) . This is a new result which was not derived previously. Let us explain the physical 
meaning of each factor restricting us to the first case. Remembering that each vortex carries an 
internal orientation of CP 1 , we can understand that the factor (Air / (g 2 c)) 2 corresponds to the 
internal orientations (see Eq. (14. 9 j) with N = 2 below) while the factor A does to the center of 
mass. Then the remaining factor A—%^- can be thought of as effective area of relative motion 
moduli of two non-Abelian vortices. Comparing this result to that of two Abelian vortices in 
Eq. (13.51) with k = 2, the extra factor 2/3 appears here, which implies that the effective area of 
non-Abelian vortices is smaller than that of Abelian vortices. Interestingly, the factor 2/3 in the 
exclusion area is different from the factor 1/2 in the Bradlow area B2 = contrary to the 

Abelian case Eq. ( 13.6ft in which the both factors coincide. We will see this in more detail around 
(Q2l> . below. 



2) k Local non-Abelian vortices with iVc = N-p = N 

Let us next consider the single (k = 1) non-Abelian vortex. In this case the moduli space is 
T 2 x CP^ -1 . The corresponding configuration has a rod and N — 1 particles trapped inside a 
rod, see Fig. H2-(ii)-(a). 

The partition function (14.41) is defined only in the region A > 47r/(g 2 c) since the number of 
eigenvalues of S(xi) is iV c and the period 27rP' should be larger than the length d of a rod: 
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5 \ 

4 
3 
2 
1 



(a) iVc = Np = 5, k = l 




(b) N t 



(i) Configurations of eigenvalues of ^(x 1 ) on slant torus. 



li- 



ft 



d 



2wR f 
(a) Nn = Np = 5, k 
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11 



d 



d 



2wR' 

(b) iV c = 1, N F = 4, fc = 2 



(ii) Configurations of rods and particles. 

Fig. 8: Examples of stant configurations on eigenvalues of ^(x 1 ) and corresponding rods and paricles. 
(i)-(a) There are 4 small kinks trapped in the region with length d. (i)-(b) There are 3 small kinks 
trapped in each interval between two regions with length d. (ii)-(a) There are 4 particles trapped inside 
a rod. (ii)-(b) There are 3 particles trapped in each interval between 2 hard rods. 



namely this bound is the second inequality in Eq. (13.21) . 

Z^r N = (^) 7V (2 7 rc)-(2^)-2^ I ^l Tyy rf 



N-l 



= W^- {dF) " A {fc) ■ <4 ' 9) 

In this simple case, we can confirm that our result agrees with an explicit integration over 
the exact metric on the non-Abelian vortex moduli space. Although the solution is not known, 
the metric with the Kahler class can be calculated [12J, [31] , to give 

ds 2 = 2ncdzdz + y^ + lftT)-^ dh i d v u .10) 

9 2 (l + |6 fe | 2 ) 2 1 ^ 

where b l are the inhomogeneous coordinates of orientational moduli CP^" 1 . Using this metric, 
we can confirm that Eq. ( 14.91) agrees with the correct partition function. 
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For general number k of the non-Abelian vortices with Ac = Np = N, we can calculate the 
partition function (14.41) up to the next leading term in the expansion in terms of 1/A 



7 N C =N ¥ =N _ (nT \kN J- 



.4 fln\ N -^ L 



(N-l)\ \g 2 cj 



k If Att 



A \ \g 2 cA 



■ (4-11) 



We can show that the coefficient Dn in Eq. (14.111) takes the form (the definition of is given 
in Appendix) 

D N (2JV-2)!! 2 8 16 128 



4tt/c/ 2 c (2JV-1)!! ' 3' 15' 35' 315' 



(N= 1,2,3, 4,5, ••■)• (4.12) 



The first (leading) term in the partition function (14. lip represents the situation that all vortices 
are separated, while the second (next leading) term implies that a pair of adjacent vortices is 
overlaped. From the partition function (14.111) . we can obtain the equation of state for dilute 
vortex gas in the thermodynamic limit A — > oo, k — > oo with fixing k/A to enough small value 



as 



PA[1-dA + o(^- \ ) =kT. (1.13) 



A \A 2 

From this we see that the second virial coefficient represents the effective area of the non- 
Abelian vortices in dilute gas. We find that the third term 0(k 2 /A 2 ) does not vanish. This 
implies that the equation of state deviates from van der Waals one, contrary to the Abelian case 
Eq. (1331). 

We find that Dn < D\ = 4^ and D N begaves as 

9 c 

i nr 47r 

°» - 5 Ht~o (414) 

for large N. We see that the non-Abelian vortices can be closer to each other than the Abelian 
vortices even though the area of individual vortex is 4r for both non-Abelian and Abelian 
vortices. Therefore we conclude that non-Abelian vortices are "softer" than Abelian vortices. 

Clearly the inequality B N < D N holds with B N = jt^- being the Bradlow area in ( 12. 51) . This 

iv g c 

inequality implies that the effective area of a vortex in dilute vortex gas is larger than the 
Bradlow area B^, which is the effective area in highly pressured gas. An intuitive understanding 
of this inequality is as follows. It is known that the internal orientations of two non-Abelian 
vortices are almost always aligned when the two vortices are approaching each other [20J (as long 
as their speed is sufficiently slow). In other words, nearby non-Abelian local vortices behave as 
if they are Abelian local vortices. 



3) k semi-local vortices with iV c = 1 and general iV F 

Finally we show that our method can be extended to the case of semi-local vortices with 
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iVp > iVc- Here we concentrate on the Nq = 1 case which indicates essential features of semi- 
local vortices. The corresponding configuration has Np — 1 particles trapped in each interval 
between k hard rods, see Fig. [SJ-(b). In this case, Eq. (14.41) reduces to 



rp \ kN F | | 

Z^7 hNF = (iz) (27tc) kNF (27tR) hNF ^jj-^— — —2ixR' (2ttR' — dk) kNF ~ 1 



2tx) v ' v k (kN F - 1) 

cT) kNF \ n J (a-^) , (4.15) 



fc(fciV F -l)! V 9 2 c 

where 1/k factor is needed since k vortices cannot be distinguished. Substituting Np = 1 into 
this result, the partition function of k ANO vortices in Eq. (13. 5p is correctly recovered. From this 
partition function, we can obtain the equation of state for the vortex gas in the thermodynamic 
limit A — > oo, k —>■ oo with k/A finite as 

P ( A -^)= kN ^ ( 416 ) 

where P is the pressure. The factor Np in the right hand side of Eq. (I4.16P appears due to the 
fact that the vortices have additional internal degrees of freedom in the case of Np > 1. 

In this case, no field theoretical result has been known yet. Especially we can show that the 
moduli space of k — 1 semilocal vortex (see Fig. [8]-(b)) is T 2 x CP^' 1 by extending the moduli 
matrix formalism (15] to T 2 . Second, as was done in [31], we can compute the moduli space 
metric by using the integration formula of the Kahler potential in [T7], to give 

9 / 4vr\ SM1 + \b k \ 2 ) - W . -. , 
ds 2 = 2vrc dzdz + 2c [A - — 1 1 ' , db l dV . 4.17 

Here b % in the second term are the inhomogeneous coordinates of CP Nf ~ x and represent moduli 
not localized around vortices. These moduli become non-normalizable on the plane C (A — > oo). 
However they contribute to the Kahler potential since we are considering the vortices on the 
compact manifold T 2 with the finite area A. From this metric we can confirm at least for k = 1 
case that Eq. (14.151) gives the correct partition function. 



5 Comments and Discussions 

We here comment on an important fact that there is a duality on the partition function (14. 4p . 
without performing the difficult integration. That is, the partition function for k < Nc is 
invariant under exchange 

{N c , g 2 ) <- (iV c , ~g 2 ) , (5.1) 

~ 47T 4tt 

with N c = k + N F -N c , — = A — . 

9 c 9 c 
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One can confirm that this duality exists between the partition function (14.91) and the one of 
the k = 1 case of (I4.15p . Also, the example of Nq = Np = 2, k = 2 is self-dual. In fact, in the 
partition function (14.81) . the latter case (A/2 < 4n/(g 2 c)) can be obtained by exchanging Am / (g 2 c) 
by A — 47r/(g 2 c) in the former case (A/2 > 4ir/(g 2 c)). These observations are very analogous 
to the duality relation between the non-commutative instantons and monopoles [32} [33] . From 
string theoretical point of view, there appears a slant torus due to the effect of the magnetic 
gauge field or equivalently the I?-field. The I?-field is a source of the non-commutativity and 
associated with the FI parameters in the effective theory on the instantons. In our case, the 
gauge coupling plays a role of the non-commutative parameter and should relate to the modulus 
of the slant torus. Better explanation on the above duality may be found in terms of the duality 
between the non-commutative vortex and domain wall in the string theoretical picture. 

Finally we would like to discuss vortices on a sphere S 2 . In this case we cannot use the 
above T-dual argument naively since the sphere consists of non-trivial U(l) fibration. However 
we here assume that the vortices on the sphere is equivalent to the vortices on the cylinder with 
finite length in the computation of the volume of the moduli space; The cylinder is topologically 
isomorphic to the sphere with two puncture points. When at least one vortex sits on a puncture 
point (the notrh or south pole) of the sphere, the dimension of the corresponding subspace is 
less than the dimension of the full moduli space. Therefore this subspace does not contribute to 
the volume of the moduli space. In contrast to the torus case, we do not need to identify the 
both sides in the x 1 -direction. Thus it is now sufficient to consider the gas of the hard rods in a 
finite segment with length 2ttR'. By using the system of the hard rods, we obtain the partition 
function of vortices in the Abelian-Higgs model on S 2 as 

Z N k ^ = 1 (cTf (2nR> - kdf (2nR) k = I (cT) k (A-^j\ (5.2) 

This completely agrees with the results in j2j H[ E]. This case also should be extendable to the 
non-Abelian and/or semi-local vortices. 

In conclusion, we have proposed a novel and simple method to compute the partition function 
of vortices at finte temperature. Our result agrees with previously known cases (13.51) and (15.21) 
of the local Abelian (ANO) vortices in the Abelian-Higgs model on T 2 and S 2 , respectively. Our 
method provides new results in more general cases of non-Abelian local vortices, ( 14.81) and ( 14. lip , 
and Abelian semi- local vortices, (14.151) . In the two cases of k — 1, Nq = Np = N and k = 1, 
iVc = 1 with general Np, we have confirmed that our results agree with explicit integrations 
over the exact moduli metrics (14.101) and (I4.17p . respectively. We have found that non-Abelian 
vortices are reduced under T-duality to soft rods with particles inside them while the Abelian 
vortices are to hard rods. 

Our results will be applied to the thermal vortex gas in the early Universe or in supercon- 
ductor. Extension to non-(or near-) critical coupling will be important to discuss more realistic 



16 



application. In particular the phase transition of vortex gas in the Abelian-Higgs model was 
discussed previously [3UE5]. Presence or absence of phase transitions in the non-Abelian and/or 
semi-local vortices is very interesting to explore. 
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A Virial Expansion 

In this appendix we give the defenition of which appears as a coefficient of the next leading 
term in the partition function ( 14. lip . We consider k non-Abelian local vortices with iVc = Np = 
N on a torus T 2 . As we have mentioned, the system can be well described from the picture of 
k soft rods with length d, each of them piercing N — 1 particles therein, on S 1 in circumference 
L = 2ttR', as shown in Fig. [5j We can carry out an integration with respect to only k parameters 
corresponding to positions of the rods and we find that the volume of the configuration space 
can be rewritten as 



Vsj. = t fd kN y= I d k(N ' 1)z 



k k\ 



max ( 1 _ 7 Yl Mn * ' 



L 

i=l 



(A.1) 



where the integration parameters z = z r ^, (r = 1, • ■ ■ , N — 1, i — 1, • • • ,k) are dimensionless and 
correspond to (N — 1) relative positions of the particles inside each rod, and their integration 
region is defined by £jv,fc = {{^r,i} |0 < < z%j, < • • • < Zn-i,% < 1, 1 < * < A:}. Here we have 
defined a dimensionless function, M^, for each rod (1 < i < k), 

M N)i = max(z lyi ,z 2 ,i- z^i+x,--- ,z N - 14 - z N - 2 ,i+iA- ZN-i,i+i), (A. 2) 

where we have defined 2V,fc+i — z r,i- All of them give the same contributions to the integration. 
Note that dimensionless parameter d/L appears in the integrand. We find kd > d^2 i M N ^ i in 
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this case of local vortices. Thus small kd/L guarantees that L > dJ^M/v^ in any point of 
the integration region and that the integrand can be expanded with respect to d/L, since the 
max-function in Eq. (1A.2j) can be ignored_|. The D ^ appears as the coefficient of the next leading 



term in this expansion of Vjv^, namely in the virial expansion of soft rods in one dimension: 

v »* - si {w^t o - b « * *<* - i)i L + ° ©) • (A - 3 » 

where we have defined Dn as an average value of Mn,i' 

Dn = ^ + ((N- l)!) 2 / dz N ~ l dz N - 1 max (z t - z[, ■ ■ ■ , z N -x - 4-1, 0) , (A.4) 

iV 70<21<-<2JV_1<1 

0<z' 1 <-<z' N _ 1 <l 

where Dn relates with -D^ by Dm = Dn^t- = D N 2nRd. We can easily carry out this integration 

g c 

in each region where a definite order of the integration parameters such as z\ < z[ < z' 2 < Z2 < 
■ ■ ■ . It is convenient to divide the integration regions into several ones whose integral become 
identical. Therefore calculation for D^ reduces to counting number of elements of the sets and 
we have succeeded this counting and obtained the result as (14.121) . 
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